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Abstract. We present the first results coming from the observation of Kepler's supernova remnant obtained with 
the EPIC instruments on board the XMM-Newton satellite. We focus on the images and radial profiles of the 
emission lines (Si K, Fe L, Fe K) and of the high energy continuum. Chiefly, the Fe L and Si K emission-line 
images are generally consistent with each other and the radial profiles show that the Si K emission extends to 
a larger radius than the Fe L emission (distinctly in the southern part of the remnant). Therefore, in contrast 
to Cas A, no inversion of the Si- and Fe-rich ejecta layers is observed in Kepler. Moreover, the Fe K emission 
peaks at a smaller radius than the Fe L emission, which imphes that the temperature increases inwards in the 
ejecta. The 4-6 keV high energy continuum map shows the same distribution as the asymmetric emission-line 
images except in the southeast where there is a strong additional emission. A two color image of the 4-6 keV and 
8-10 keV high energy continuum illustrates that the hardness variations of the continuum are weak all along the 
remnant except in a few knots. The asymmetry in the Fe K emission-line is not associated with any asymmetry 
in the Fe K equivalent width map. The Si K maps lead to the same conclusions. Hence, abundance variations do 
not cause the north-south brightness asymmetry. The strong emission in the north may be due to overdensities in 
the circumstellar medium. In the southeastern region of the remnant, the lines have a very low equivalent width 
and the X-ray emission is largely nonthermal. 
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1. Introduction 

The remnants of recent stellar explosions give out a large 
part of their luminosity in X-rays. The interaction of the 
high velocity supernova material (ejecta) with the ambient 
medium gives rise to a forward shock moving into the am- 
bient medium and a reverse shock propagating back into 
the ejecta IjMcKee 1974|l . Both shocks compress and heat 
the matter. Thus the X-ray emission arises from two adja- 
cent media, different in composition, density and tempera- 
ture: the shocked ejecta and the shocked ambient medium. 

At an age of nearly 400 years, Kepler is a young super- 
nova remnant (SNR) whose X-ray emission is still dom- 
inated by a thermal spectrum characteristic of shocked 
ejecta. It is about 200" in angular size. At a distance of 
about 4.8 ± 1.4 kpc ( |Reynoso fc Goss 1999| l, this corre- 
sponds to an angular radius of about 2.3 ± 0.7 pc. This 
shell-type remnant shows a strongly asymmetric bright- 
ness distribution in X-rays. Well correlated optical and 
infrared observations also exhibit an asymmetric emission, 
with the remnant being notably brighter in the north and 
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northwest with some structures near the center of the rem- 
nant (Bandiera & van den Bergh 1991; Douvion et al. 
2001). The optical knots in Kepler's SNR, which have a 
low expansion velocity (Bandiera & van den Bergh 1991), 
a high density and nitrogen overabundances (Dennefeld 
1982; Blair et al. 1991), are attributed to a circumstellar 
origin. The measurement in a Balmer-dominated knot of 
a large [Nil] line width favors a nonthermal broadening 
in a cosmic-ray precursor (SoUerman et al. 2003). As for 
the radio emission, it has the same shell-like morphology 
as the X-rays with the same north-south brightness asym- 
metry (Matsui et al. 1984; Dickel et al. 1988). Otherwise, 
Kepler's SNR's type is not clearly established. On one 
hand, its position above the Galactic plane (~ 600 pc) 
and its geometrical similarities with Tycho's SNR favor 
a type la progenitor l|Smith et al. 19891 whereas, on the 
other hand, optical and infrared observations which re- 
veal the presence of circumstellar material suggest a type 
II progenitor. The light-curve does not allow to distin- 
guish between a type I IjBaade 1943(1 and a type II-L 
| |Doggett fc Branch 1985| ) classification. 

Previous X-ray studies have been made on images 
and spectra, but always separately. The best images 



01 Kepler s biNK were obtained with the hmstem HKi 
(Whi te fc Long 19 83; Matsui et al. 1984) and ROSAT 
HRI ( [Hughes 1999| |. The upper Umit of the energy range 
did not exceed 4 keV and no spectro-imaging was avail- 
able. Spectral studies were more numerous. The SSS on 
board the Einstein satellite revealed line emission from the 
He-like species of Si, S and Ar HBecker et al. 1980|l . The 1 
to 10 keV spectra obtained with the EXOSAT and Ginga 
observatories have shown a strong iron Ka emission line at 
nearly 6.5 keV ll Smith et al. 19891 iHateukade et al.TMiji . 
ASCA data confirmed these facts with a better energy res- 
olution | |Kinugasa fc Tsunemi 1999| ). From these spectra, 
several models have been tested coupling non-equilibrium 
conditions with the hydrodynamic evolution of the SNR 
(Hughes fc Helfand 1985; Decourchelle & Ballet 1994; 
Rothenflug et al. 1994; Decourche lle et al. 2000). 

A more specific model ( Borkowski et al. 1994|l was ap- 
plied following Bandiera's idea l|Bandiera 1987|l in which 
the progenitor was a massive runaway mass-losing star 
ejected from the Galactic plane. A specific review on 
Kepler's SNR can be found in Decourchelle (2000). 

The new generation satellites XMM-Newton and 
Chandra mark a new era in observing capabilities. Indeed, 
the onboard technology is made of CCD detectors as- 
sociated with telescopes allowing combined good spec- 
tral and spatial resolution with in addition the capabil- 
ity to collect a large number of photons. Furthermore, 
both satellites are equipped with slitless gratings giving 
access to high resolution spectroscopy for point or slightly 
extended sources particularly around 1 keV (around the 
Fe L complex). These improvements enable us to analyse 
thoroughly the spectra of small-scale structures of a SNR 
as well as to study maps in a given energy band. This kind 
of work has been done for several SNRs: Cas A (Hughes et 
al. 2000; Hwang et al. 2000; Bleeker et al. 2001; Gotthelf 
et al. 2001; Willingale et al. 2002, WilHngale et al. 2003), 
Tycho (Decourchelle et al. 2001; Hwang et al. 2002), the 
Crab Nebula (WilHngale et al. 2001). 

Thanks to the EPIC imaging spectrometers of XMM- 
Newton observatory, it is henceforth possible to carry on 
studying more precisely Kepler's SNR. In this paper, we 
intend to draw up maps of the brightest X-ray emission 
lines and to lay emphasis on the high energy continuum. 
Together with this global approach, we present early re- 
sults coming from a local spectral study of the forward 
shock. 

2. Data reduction 

Kepler's SNR was observed on 2001 March 10 using the 
XMM-Newton EPIC and RGS instruments. The result of 
the RGS analysis is complicated by the large spatial extent 
of the source and will be the purpose of a forthcoming 
paper. Here, we focus only on the photons collected by 
the EPIC instruments : MOSl, M0S2 and pn. The thick 
filter has been used to limit the number of low-energy 
photons and hence pile-up. For the pn camera, we have 
chosen the large window mode in order to limit the out- 
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Table 1. Flux loss and pile- up rate in the worst case cor- 
responding to the 0.0086 cts/frame/pixel of MOS. 



Si_K SK Fe K 

MOSl 1.8606 ± 0.0004 2.4465 ± 0.0010 6.4655 ± 0.0050 

MOS2 1.8540 ± 0.0004 2.4374 ± 0.0010 6.4374 ± 0.0050 

pn 1.8471 ± 0.0004 2.4304 ± 0.0009 6.4481 ± 0.0028 

ASCA 1.854 ± 0.002 2.445 ± 0.003 6.46 ±0.01 

Table 2. Fitted centroid energies (in keV) of three promi- 
nent line blends in Kepler's SNR after CTI correction as 
from the integrated spectrum. The errors are in the range 
Ax^ < 2.7 (90% confidence level). 



of-time events. The obtained data quality is very good 
(few flares). We have used the Science Analysis System 
(SAS, version 5.3) for data reduction. 

We have discovered in the M0S2 data the pres- 
ence of a few bright pixels in the central CCD which 
shift and broaden the spectrum of the whole column 
(in the readout-out direction) containing these bad pix- 
els. They existed before the launch of XMM-Newton but 
were never detected. Consequently, we have removed the 
entire columns RAWX=273, 276, 277 and the column 
RAWX=383 from RAWY=311 to RAWY=600 in the 
MOS 2 events list. 

Because the source is bright, it is necessary to esti- 
mate pile- up HBallet 1999|l . The maximum brightness (in 
cts/frame/pixel) is equal to 0.0086 for MOS and 0.0048 
for pn. We can estimate the flux loss and the pile-up rate 
at the maximum brightness^. The results are shown in 
Tabled They show that pile- up is not enough to alter 
our results, although the true images would show slightly 
more contrast. 

To create light curves, spectra and images, we have 
selected the single, double, triple and quadruple events 
(pattern<12) for the MOS cameras, and the singles and 
doubles (pattern<4) for the pn camera (with FLAG=0 for 
all cameras). The data were cleaned from the few flares 
using the light curves between 10-12 keV for MOS and 
12-14 keV for pn. Then, the exposure time of the obser- 
vation amounts to 30 ks for each MOS camera and to 27 
ks for the pn camera. Although MOS data were corrected 
for Charge- Transfer Inefficiency (CTI), a spectral shift de- 
pending on energy remains (see Table ^ . The fitted cen- 
troid energies from the EPIC instruments and the ASCA 
data (Kinugasa fc Tsunemi 1999) emphasize that a single 
additive or multiplicative factor for each data set cannot 
be found so that the individual lines (e.g.. Si K, S K and 

^ The result depends on the proportion of X-rays creating a 
charge pattern with i pixels called a,; in Ballet (1999). Here 
ai = [0.708; 0.233; 0.021; 0.038] (from the data, averaged over 
energy). 
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Fig. 1. The EPIC spectra obtained for the entire remnant 
(MOS below and pn above). The energy bands in which 
we are interested are indicated in Table 01 

Fe K) can be made mutually consistent for the various 
instruments. Moreover, since there are no absolute posi- 
tion for the line blends due to non-equilibrium effects, we 
cannot re-calibrate these line blends at a reference energy. 

Kepler's SNR is contained within the central CCD of 
the MOS cameras and spread over 4 CCDs of the pn cam- 
era resulting in gaps between each CCD. Below 4 keV, 
we have only used the MOS data since we have enough 
statistics to obtain precise images without gaps. Above 4 
keV, the MOS data have been combined with the pn's to 
remedy the lack of statistics at higher energy but intro- 
ducing gaps in the images. Note that the spectral analysis 
has been carried out over 0.3-10 kcV and 0.3-12 keV for 
MOS and pn, respectively. 

In this paper, we did not use a local background to be 
subtracted (in the spectra as well as in the images) since 
Kepler's SNR is so bright that the background is con- 
taminated by the wings of the PSF over the entire field 
of view (FOV). Thus, we have only subtracted a 400 ks 
astrophysical and instrumental background (D. Lumb et 
al. 2002) which has been renormalized in the 10-12 keV 
band for MOS and 12-14 keV band for pn over the full 
FOV. This background has been obtained using a thin 
filter which is different from the thick filter of our ob- 
servation. However, it does not affect the lower energy 
spectrum because it is more than one order of magnitude 
below Kepler's SNR spectrum at 0.3 keV. For the follow- 
ing spectral studies (see SectE)), we have used the method 
which consists of calculating a weight for each photon to 
correct for the vignetting effect. The chosen on-axis re- 
sponse files are m*_tliickv9q20t5r6^11_15.rsp for MOS 
and epn_lw20_sdY9_thick.rsp for pn. 

The EPIC spectra extracted from the entire remnant 
are shown in Figure^] These X-ray spectra are expected to 
be the combination of at least two optically thin thermal 
continua plus emission lines coming from highly ionized 
atoms of the ejecta. The first continuum arises from the 
shocked ejecta and the second continuum comes from the 
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Table 3. Energy cuts for images. 



shocked ambient medium. A nonthermal component (due 
to synchrotron process or nonthermal bremsstrahlung) 
may also be present. Chiefly, He a features of silicon, 
sulphur, argon, calcium and the iron Ka line dominate 
the EPIC spectra from about 1.4 keV. Below, numerous 
transitions of the iron L lines are blended with magne- 
sium, neon and oxygen K lines. At first sight, the over- 
all spectrum resembles Tycho's (Decourchelle et al. 2001). 
However, it is noteworthy that the Fe L complex and the 
Fe K line respectively dominate over the Si and the Ca 
lines, unlike in Tycho. 

The emission-line images are formed from these spec- 
tra by selecting narrow energy bands. Table 13 gives the 
energy cuts chosen for these images. The width of the 
Si K and Fe K energy bands have been optimized (best 
signal-to-noise ratio) by a method explained in Appendix 
lA.ll We also cleaned the images contaminated by con- 
tinuum contribution with an interpolation based method 
which is described in Appendix IA.2I Above all, this was 
done for the Fe K image using the surrounding continuum, 
which is well represented with an exponential model. The 
method was also applied to the Si K using a power-law 
model. However, one should keep in mind that this method 
is approximate since the base level around the Si K line 
does not correspond exactly to the continuum level. This 
method was not applied to the Fe L line owing to the fact 
that the continuum level is not well known around this 
line. 

The mapped lines are those that have the highest 
signal-to-noise ratio above the continuum: Si Ka, Fe L 
and Fe Ka. We do not show maps of fainter lines of S, Ca, 
and Ar that are present in the integrated spectrum. They 
are roughly similar to the map of the Si K line emission. 

3. Emission-line images 

Figure |21 (top) shows Fe L and Si K adaptively smoothed 
images. All adaptive smoothings were made with the as- 
mooth package of the XMM SAS. FigureEl (bottom) shows 
their respective deconvolved images to better illustrate the 
main features of the SNR. Deconvolved images have been 
constructed using the MOS PSF at 1500 eV (estabhshed 
from the calibration files) and the method is the regu- 
larized gradient method with 40 iterations, assuming the 
noise to be Poisson l|Starck et al. 1998|l . The FWHM of 
the PSF is 4.4 arcseconds at 1500 eV. 
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Fig. 2. Top-left panel: MOS vignetting corrected count rate image of Fe L (780-860 eV band). The data are smoothed 
with the same template used for the Si K image (Top-right). Bottom-left panel: Fe L deconvolved image in the 
same band. Top-right panel: MOS vignetting corrected count rate image of Si K (1755-1940 eV band). The data are 
adaptively smoothed to a signal-to-noise ratio of 5. Bottom-right panel: Si K deconvolved image in the same band. 
All images use a square-root scaling. The small black circle in both deconvolved images indicates the center chosen to 
construct the radial profiles. 



In young SNRs, we expect to see at least two differ- 
ent emission zones coming on one hand from the shocked 
ejecta and on the other hand from the shocked ambient 
medium. The X-ray emission, crudely proportional to the 
density squared, mainly arises from the shocked ejecta 
which are usually denser than the ambient medium and 
enriched in heavy elements. Thus, the brightest emission 
lines give information on the ejected supernova matter. 

The comparison between the spatial distribution of Fe 
L and Si K indicates an overall morphological agreement 
both in the almost spherical diffuse emission and in some 
of the brightest regions. Indeed, the north-south bright- 



ness asymmetry of the ring, the bar, the elbow, the west- 
ern ear, and the eastern eye are shared by both images. A 
closer look reveals that some differences exist notably at 
the eastern ear where there seems to be more Si K than 
Fe L and at the hand where it is the contrary. We have 
constructed separately the radial profiles of the northern 
(10°-170°, from west anti-clockwise) and southern (175°- 
360°) parts of Kepler's SNR from the images. These radial 
profiles (Fig. emphasize the strong contrast between 
the north (left panel) and the south (right panel). The 
north emission is exclusively dominated by the northern 
ring whereas at the south, a large part of the inner emis- 
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Fig. 3. Left panel: MOS+pn radial profiles of Fe L (solid line), Si K (dashed), Fe K (dotted) and 4-6 keV continuum 
(dash dot) between 10° and 170° from west (northern ring). The coordinates of the center (17h30m41s,-21°29'23") 
are taken from Hughes (1999). It is marked by the small black circle on Figure El The radial profiles were built from 
the vignetting corrected count rate images of Fe L, Si K and 4-6 keV continuum (background subtracted) and Fe K 
(continuum and background subtracted) with a step radius of 4.4". Each radial profile has been renormalized to its 
maximum. Right panel: same as left panel but between 175° and 360°. 



sion is due to the eastern eye. We also see that both Fe 
L and Si K spread out further in the south than in the 
north. These radial profiles show that the Si K emission 
peaks at a larger radius than does Fe L: slightly in the 
north and rather more in the south. It shows that there 
is no inversion of the Si and Fe layers as detected in Cas 
A (Hughes et al. 2000; Hwang et al. 2000; WiUingale et 
al. 2002). These clues agree with an onion layer structure 
(at least for iron and silicon) which would have persisted 
after the explosion. 

The comparison between the Fe L and Fe K emission 
can give information on the temperature structure inside 
the remnant and on the ionization timescale. Figures |31 
and 01 illustrate that the Fe K emission peaks at a smaller 
radius than the Fe L emission in the north as well as 
in the south. This tells us that the inner ejecta should 
be hotter than the outer ejecta, but not directly by how 
much. Kinugasa & Tsunemi (1999) have modelled the in- 
tegrated spectrum of Kepler's SNR with the ASCA data 
using two non-equilibrium thermal components with the 
same abundances for both components but with different 
temperatures. They found a temperature of 0.36 ± 0.04 
keV and 2.26 ± 0.15 keV for the Fe L and Fe K, respec- 
tively. Even though the Fe K was underproduced with this 
two-component model, this provides evidence for a strong 
temperature gradient (a lower limit temperature ratio of 
~ 6) in the shocked ejecta medium turned inwards of the 
remnant. We note that an ionization effect would work in 
the opposite direction. The ratio between the Fe K and 
the Fe L line strengths should increase when going away 
from the reverse shock as Fe gets more ionized. 

The temperature increase inwards must be the result 
of the remnant's hydrodynamics. Hydrodynamical mod- 
els give the profiles of pressure P and mass density g as a 
function of the radius. We note that the profile of P/ g is 
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Fig. 4. Fe L contours of the deconvolved image superim- 
posed with the MOS-|-pn Fe K vignetting corrected count 
rate image (continuum and background subtracted) in the 
6180-6169 eV band. The data are adaptively smoothed to 
a signal-to- noise ratio of 5. The intensity scaling is square- 
root. 

equivalent to the temperature profile if the mean molecu- 
lar weight remains constant, i.e. if the abundance of the 
dominant elements does not vary in the remnant. In the 
following paragraphs, we assume constant elemental com- 
position but we should keep in mind that any variation 
of the ejecta composition of the dominant species should 
alter our interpretation. 

Models of SNRs evolution depend on the initial outer 
density profiles in the ejecta and the ambient medium. We 



tirst assume that the supernova exploded into a constant 
density medium. Two simple cases have been under con- 
sideration in the framework of spherical models. On one 
hand, the ejecta density profile may follow a power-law 
distribution with a constant core I) Chevalier 1982|l which 
is appropriate to a type II classification IjArnett 1988|l . 
In this hypothesis, contrary to the observations, the 
ejecta are hotter near the contact discontinuity (inter- 
face between the shocked ejecta and the shocked ambient 
medium) than at the reverse shock. On the other hand, 
an exponential density profile aims to represent a type 
la progenitor IjDwarkadas fc Chevalier 1998|l . In this lat- 
ter hypothesis, the shocked ejecta are gradually colder to- 
wards the contact discontinuity as is observed, but the 
temperature profile is much too flat to account for the 
required two temperatures of the Fe K and Fe L ejecta 
(temperature ratio of ~ 6). The temperature gradient of 
the reverse shocked region can be larger when the expo- 
nential ejecta density structure is running into a circum- 
stellar medium ( Dwarkadas fc Chevalier 1998|l and is then 
compatible with the observations^. 

In the remnant of Tycho's type la supernova, the 
Fe K emission peaks at a smaller radius than the Fe L 
(Decourchelle et al. 2001; Hwang & Gotthelf 1997) as is 
observed in Kepler's SNR. Both remnants share the prop- 
erty of a higher temperature toward the center in the 
shocked ejecta medium with a temperature ratio of at least 
2 (Hwang et al. 1998; Kinugasa & Tsunemi 1999). In the 
framework of an exponential ejecta profile, this requires 
an interaction with a circumstellar medium. There is no 
observational indication for such a circumstellar medium 
around type la supernovae but as discussed in Dwarkadas 
& Chevalier (1998), it has been proposed that in some bi- 
nary models of type la the mass loss of the companion can 
modify the environment of the progenitor. 

In Kepler's SNR, there is evidence for the presence 
of a circumstellar material which has been interpreted in 
terms of a massive progenitor (Bandiera 1987; Borkowski 
et al. 1992). If the supernova was a type lb as suggested 
by Borkowski et al. (1992) and has a power-law outer den- 
sity profile, the observed higher temperature inwards the 
remnant would require that the reverse shock has already 
entered the central plateau. This condition requires that 
the ejecta mass Mej is lower than 3Mq taking the pa- 
rameters^ derived from the best fit to the ASCA+XTE 
data (Decourchelle et al. 2000). If Mej > 3Mq, the re- 
verse shock would still be in the power-law envelope, and 
would give the wrong behavior in Fe K/Fe L. 



^ We note that a power-law ejecta density structure running 
into a circumstellar medium would also lead to a temperature 
slope compatible with the observations but with a too low tem- 
perature to account for the Fe K line (Decourchelle & Ballet 
1994). So this case is rejected and no more considered in the 
following. 

^ EsNiMQ/Mcif/^ = 0.53 10^^ ergs and po = 0.42 amu 
cm~"^ for a power-law index of the initial outer density profile 
in the ejecta n — 9 and a flat ambient medium. 
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Fig. 5. MOS-|-pn vignetting corrected count rate im- 
age of the 4-6 keV continuum. The data are adaptively 
smoothed to a signal-to-noise ratio of 5. The intensity 
scaling is square-root. The spectrum of the region Al 
is shown in Figure |H1 The center of the region Al is 
(17h30m47s,-21°30'39"), its radius is - 14" and it is 
^ 105" away from the SNR center. 

The condition that the reverse shock has entered in 
the plateau sets an upper limit on the expansion rate of 
the SNR. Indeed from the self-similar model of Chevalier 
(1982), the expansion parameter is m = (n — 3)/n oc 
log i?/ log i assuming a power-law density profile of in- 
dex n expanding into a flat ambient medium. At most 
rt=12 imposes the upper limit m = 0.75 for a flat ambient 
medium. This value is contradictory with the results of 
Hughes (1999) who measures a dynamical evolution near 
free expansion (m ~ 0.93). However, it is consistent with 
the highest value m ^ 0.65 found from the radio measure- 
ment (Dickel et al. 1988). 

4. High energy continuum 

As we see in Figure |S1 the 4-6 keV continuum emission 
globally looks like the emission-line images except in the 
southeast where there is a clear limb-brightened profile 
which broadly resembles the 21 cm radio image (Dickel 
et al. 1988; DeLaney et al. 2002). Radial profiles (Fig.EJ 
and the three color image (Fig. left panel) indicate that 
the continuum emission occurs at a larger radius than all 
the line emission, especially in the southeast. We surmise 
that the outermost continuum emission is the signature of 
the blast wave. The north emission is the strongest. The 
southeastern side (near the eastern ear) defines a ridge of 
emission dominated by the presence of a few bright fea- 
tures at XMM-Newton spatial resolution. These features 
seen with XMM-Newton correspond to thin curved shocks 
observed in the Chandra image. The spectrum of the re- 
gion Al (cf. small circle in Fig.^J will be analysed in Sect. 




Fig. 6. Left panel: Three color image of the Si (blue), Fe K 
(red), 4-6 keV high energy continuum (green). All images 
are vignetting corrected and are smoothed with the tem- 
plate created from the adaptively smoothed count image 
of the Fe K -I- 4-6 keV (signal-to- noise of 5). Right panel: 
MOS-|-pn two color image of the high energy continuum: 
4-6 keV map in green, 8-10 keV map in red, common ar- 
eas are in yellow. Both count rate images are background 
subtracted, vignetting corrected and are smoothed with 
the template created from the adaptively smoothed 8-10 
keV band count image (signal-to- noise of 5). Maps use a 
square-root scaling with minimum set to 30% of maxi- 
mum. 

IHI A diffuse emission region exists at the center of Kepler's 
SNR toward the eastern eye. The southwestern edge is less 
pronounced than the other outer rims. 

Figure ini (right panel) shows that the 4-6 keV and 8-10 
keV bands correspond well and that the continuum emis- 
sion essentially arises from an outer ring. It also illustrates 
how the hardness of the continuum changes all along the 
remnant and gives locally an indication of the shape of 
the high energy spectrum. Regions where the continuum 
is flat (in orange) are found in the north, at the shoulder 
and in the southeast. It is possible that these regions are 
responsible for the higher energy emission (up to 20 keV) 
discovered by RXTE UDecourchelle fc Petre 1999|l . 

The northern part of Kepler's SNR is the brightest 
in the silicon and iron maps as well as in the high en- 
ergy continuum map (see for example Fig.lHl left panel). 
At first sight, it may be due to a density enhancement in 
the ambient medium toward the north. Since brightness is 
proportional to the square of the density, such an enhance- 
ment would explain the strong brightness of the northern 
rim in all the different energy bands that we have selected 
to construct our images. The southeastern emission (from 
the shoulder to the southern ring), which is faint in the 
emission lines, may essentially be due to a low ambient 
medium density. To check these assumptions, it is possi- 
ble to map the ratio of the line emission to the underlying 
continuum, i.e. to construct equivalent width maps. 



5. Equivalent width images 

We focus on the Fe K and Si K equivalent width (here- 
after EQW) for which it is possible to define an adjacent 
continuum. The EQW map reveals the ejecta abundances 



but depends also on the temperature and on the ioniza- 
tion age. Note that it also assumes that the underlying 
continuum is related to the lines. To construct an EQW 
image, we have divided the image in the line (background 
and continuum subtracted) by the underlying continuum 
(background subtracted) [see Appendix IA.3j . 

Figure (left panel) shows the Fe K EQW. What 
we see is the contrast between the center of the remnant 
and the exterior. The former is associated with the ejecta 
(high EQW) and the latter with the forward shock (weak 
EQW). Hence, this confirms that the iron originates from 
the ejecta. We note that the Fe K ejecta are mainly found 
at the shoulder, the eastern eye and below the hand (high 
EQW). On the other hand, the southeastern region is ex- 
tremely faint in Fc K (weak EQW) but bright in the high 
energy continuum. 

When focusing our attention in the inner part of the 
ejecta (in red. Fig. [7| left panel), we see that there is no 
obvious correlation between the EQW map and the flux 
in the Fe K ejecta. In particular the northern rim is no 
longer dominant in brightness. This is an indication that 
abundance variations are not the cause of the north-south 
brightness asymmetry and implies some differences in the 
density of the ambient medium around Kepler's SNR. This 
point is consistent with the presence of dense optical knots 
in the northern ring (Bandiera fc van den Bergh iMTl ). 

Figure [71 (right panel) shows the Si K EQW. The con- 
trast between the center of the remnant and the exte- 
rior is less pronounced than in the Fe K EQW. The Si 
K ejecta are found towards the shoulder, near the south- 
western edge, at the crossing ring and at the western ear 
(in yellow and red). 

Firstly, it shows afresh that the Si K ejecta have been 
dispersed to a higher radius than the Fe K ejecta, notably 
in the south. Secondly, the northern rim is not strong and 
confirms again that the north-south brightness asymmetry 
is not caused by abundance variations. As for the south- 
eastern side of the remnant, it presents a lower contrast 
than in the Fe K EQW because the Si K extends to a 
larger radius than the Fe K as shown in the radial profiles 
(Fig. 131 right panel). While the Si K ejecta emission is 
present in the western ear, the eastern ear does not show 
evidence for the presence of ejecta. 

The search for prominent interstellar Hi features 
around Kepler's SNR has lead to negative result (Reynoso 
& Goss 1999). Thus, it is likely that the presence of over- 
densities in the north is circumstellar in origin. This is 
supported by the detection of slow-moving optical knots 
{r-^ 200 km s~^) there that have nitrogen overabundances 
higher than the solar value by a factor 4 (Dennefeld 1982; 
Blair et al. 1991). 

6. Spectral study of the forward shock 

In this section, we intend to characterize the spectrum 
of the southern regions dominated by strong high energy 
continuum (4-6 keV and 8-10 keV) and weak EQW of 
Si K and Fe K. For these regions, we attribute the X- 
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Fig. 7. Left panel: MOS+pn equivalent width image of the Fe K emission. This image uses a square-root scahng. 
Colorbar is in keV: purple-blue=0-l, green=l-2.7, yellow-orangc=2.7-4.5, rcd=:4.5-6.7, white=6.7-7.4. The ratio refers 
to images smoothed with the template created from the adaptively smoothed count image of the Fe K + 4-6 keV 
(signal-to-noise of 19). The superimposed contours (sqrt scale) come from the MOS Fe L deconvolved image. Right 
panel: MOS equivalent width image of the Si K emission. This image uses a square- root scaling. Colorbar is in keV: 
blue=0-0.2, green=0.2-0.7, yellow-orange=0. 7-1.1, red=l. 1-1.7, white=l. 7-2.0 The ratio refers to images smoothed 
with the template created from the adaptively smoothed count image of the Si K -I- 1.42-1.50 keV (signal-to-noise of 
19). The superimposed contours (sqrt scale) come from the MOS Si K deconvolved image. Both images indicate the 
location of the bright region Al seen in the 4-6 keV high energy contimmm image (see Fig. El- 
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Fig. 8. EPIC spectra of the southeastern part of Kepler's SNR (MOSl in green, M0S2 in red and pn in black). Left 
panel: Best- fit obtained with the single NEI model where all the elemental abundances are fixed to solar abundances 
(see Sect.lOand TableU first model). Right panel: Best-fit obtained with the SRCUT+NEI model (see Sect. 16.2.21 
and Table |S1 second model) . For this model, we have illustrated the decomposition of the different components for 
MOSl only. Both images shows the EPIC spectra rebinned to 5cr to better distinguish the different components but 
all the spectral analysis was made with a threshold at 3(t. 



ray emission to the shocked ambient medium located just 
behind the blast wave. We will not focus in detail on the 
modelling of the contribution of the shocked ejecta which 
arises in our spectrum. Note that a detailed study of the 



forward shock all along the southeast rim would require 
Chandra spatial resolution. 

Here, we will present the first attempt to character- 
ize the spectrum associated with the forward shock in 
Kepler's SNR. For this purpose, we have extracted the 
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Table 4. Best-fit parameter values for the different single 
non-equilibrium models using the MOS and pn data all 
at once. The first model refers to the NEI model with 
all elements fixed to solar abundances. The second model 
refers to the NEI model with all the elemental abundances 
(C, N, O, Ne, Mg, Si, S, Ca, Fe, Ni) freely fitted excepted 
H and He which were frozen to solar abundances. The 
second model leads to abundances of C, N, O, Ne, Ca 
near zero. The errors are in the range Ax^ < 2.7 (90% 
confidence level) on one parameter. Quantities in square 
brackets are held fixed at the given value. 
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Si (Abo) 
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0.8 (0.4-1.2) 
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0.11 (0.09-0.13) 
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0.35 (0.33-0.37) 
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Table 5. Best-fit parameter values for the different two 
components models using the MOS and pn data all at 
once. The subscript 1 refers to the hot component (i.e. 
shocked ambient medium) and the subscript 2 to the cold 
component (i.e. shocked ejecta). The errors are in the 
range Ax^ < 2.7 (90% confidence level) on one param- 
eter. Quantities in square brackets are held fixed at the 
given value. 



EPIC spectra (Fig.|Hll of the bright region Al in the south- 
eastern area (see Fig.O. These spectra have a continuum 
extending till 10 keV with little line emission (Fig. ISJ'*. 
Though spatial redistribution of emission from the bright- 
est regions in the north could contaminate the region Al, 
we find that regions one arcminute from the center of 
the region Al contribute only ~ 5% of the flux in Al 
at 1.5 keV. Moreover, the spectra of other regions along 
the southeastern rim are qualitatively similar to Al's. 

6.1. Single non-equilibrium model 

Firstly, we have supposed that the X-ray emission of 
the region Al came entirely from the shocked ambient 
medium. We applied a single non-equilibrium (NEI) model 
to the region Al within XSPEC (version 11.1, Arnaud 
1996) freezing all the elemental abundances to solar abun- 
dances (given by Anders & Grevesse 1989). In addition to 
the elemental abundances, this NEI model embraces the 
following parameters: the temperature kT and the ioniza- 
tion age T. We define r as J^^ riedt, where Ue is the electron 
density and where tg and to are respectively the age of the 
SNR and the time when the gas was shock-heated by the 
blast wave (setting zero to the explosion date). 

With all the elemental abundances held fixed at so- 
lar abundances, absorption along the line-of-sight, tem- 
perature and ionization age were freely fitted. This results 
in a poor fit {x^— 1287 for 607 dof) which overproduces 
the lines from intermediate-mass elements (O, Ne, Mg) 
and underproduces the Si line as shown in Figure |H1 (left 
panel). This shows that at least the Si comes from the 
ejecta. Figure |H1 (left panel) shows also that the centroid 
energy position of the Si and S lines differs when compar- 

* We grouped the channels of these spectra until we have a 
signal-to-noise ratio greater than Scr in each bin. 



ing the model and the data. The best-fit parameters are 
given in Table 21(Solar abundances model). 

Then, we thawed all the elemental abundances (C, N, 
O, Ne, Mg, Si, S, Ca, Fe, Ni) except H and He which were 
frozen to solar abundances. The fit (not shown) is much 
improved (x^= 693 for 596 dof) provided that the abun- 
dances of the elements C, N, O, Ne, Ca are close to zero. 
The best-fit parameters are given in Table 0] (Free abun- 
dances model). The Si line, better represented, requires a 
slightly larger value than solar abundance while the S and 
the Fe are required to be subsolar. It seems that the Fe is 
necessary for representing some portion of the low energy 
spectrum, although we do not have the spectral resolu- 
tion there with the EPIC instruments to clearly identify 
individual Fe L lines. 

Even though the values of the temperature and the 
ionization age are reasonable, such a model is not con- 
sistent with the composition of the ambient medium. 
Nevertheless, it shows that the spectrum of the region 
Al can be reproduced formally only with a continuum 
and lines of Si, S and Fe. So, we replaced the continuum 
arising from the H and He by a second component. 

6.2. Two components models 

We introduce a second component to allow for the dif- 
ferent conditions for the lines and for the continuum. On 
one hand, the lines are attributed to the ejecta and thus 
we introduce a plasma composed only of Si, S and Fe at 
low temperature. This constitutes the cold component for 
which we fixed the Si abundance to solar abundance and 
we only thawed the S and Fe elements while fixing all the 
others including H and He to null solar abundances for 
the cold component (henceforth, what counts is only the 



abundance ratio with regard to bij. the temperature kl2 
and the ionization age T2 were freely fitted. On the other 
hand, the continuum is attributed to the shocked ambient 
medium and is associated to a hot temperature plasma. 
For this hot component, we consider two possibilities: an- 
other NEI model or an additive nonthermal component. 

6.2.1. Thermal model 

In this model, the shocked ejecta are modelled by a first 
NEI at low temperature and the shocked ambient medium 
by a second NEI at high temperature. All the elemental 
abundances were fixed to solar values. The temperature 
kTi and the ionization age ri were freely fitted. 

The best- fit parameters are given in Tableland show 
that the temperature is high (~ 2.7 keV) and the ioniza- 
tion age very low (~ 10* s cm^"^) for the shocked ambient 
medium. This is very similar to what was found in Tycho's 
SNR (Hwang et al. 2002). Indeed, in Tycho's SNR the tem- 
peratures of the shocked ambient medium (corresponding 
to Behind Rim 1,2,3,4 in Hwang et al. 2002) are around 
2 keV and the ionization age has an upper limit in the 
range of a few 10* s cm^^. 

The parameters derived from our spectral fit give con- 
straints on the density via the emission measure and the 
ionization age. In the following, we check the consistency 
between these two estimates. 

A mean ionization age f can be calculated from sim- 
ple hypotheses made on our extraction region (region Al). 
For that region, we assume that the integrated volume 
along the line-of-sight is a cylinder completely filled with 
shocked ejecta and shocked ambient medium. We define 
the mean ionization age t = ^hetg where ^ is a parameter 
depending on the spatial distribution, he is the mean elec- 
tron density which depends on the filling factor e, i.e. the 
volume fraction occupied by each medium (shocked ejecta 
and shocked ambient medium) in the region Al, and ts is 
the age of the SNR. Hence, f is a function of e and ^. In 
Appendix^ we described in detail how to compute f and 
we give general formulae that can be applied to any SNR. 
Here, we just give a numerical application of the formula 
(|B.6|I for the region Al in Kepler's SNR. 

For the shocked ambient medium, we derive the fol- 
lowing mean ionization age: 

fi = 3.7^^6 10'" ^ e^'^" s cm-3 (1) 
and for the shocked ejecta: 

h - 1.5tEJ j 10' e s-'^' s cm-3 (2) 

Now, we shall compare the ionization age given by 
Tableland the mean ionization age f in our region both 
for the shocked ambient medium and the shocked ejecta. 
The equality between these two estimates must constrain 
the value of ^ e~^/^, i.e. essentially the value of the filling 
factor e. 

For the shocked ambient medium, the lower limit on 
fi is 3.1 10^° ^ e^^^^ s cm-3 (cf. Eq.^. Even if ei = 1 



(i.e. the shocked ambient medium mis all the volume 01 
the projected region of Al), we see that the mean ioniza- 
tion age fi cannot be lower than 3.1 10^" ^ s cm^'^. Since 
f ~ 1/3 is a reasonable choice (cf. Table IBAI in Appendix 
Ib)| , we see that this lower limit value is inconsistent by two 
orders of magnitude with the value ri = 10* s cm~'^ given 
in Table m This demonstrates that the model represent- 
ing the shocked ambient medium is not consistent. If the 
shocked ambient medium does not fill the entire volume 
of the projected region (Case 2 in Appendix IB.2|I . it in- 
creases the lower limit on fi which is already inconsistent. 

Yet, the Chandra image of the 4-6 keV continuum 
shows that the bulk of the X-ray emission coming from 
the region Al arises from filaments corresponding to thin 
shocks which extend all along the southwest. Then, from 
a morphological view point, it seems that this X-ray emis- 
sion arises from a recently shocked ambient medium con- 
centrated in a small fraction of the entire volume of the 
projected region. In that case (cf. Appendix IB.3p . Eqs. 
(jRTIl and (IrU leads to n = l.etj;? lO^o e^^ s cm'^. 
The equality between this mean ionization age and ri (cf. 
Table 0) implies that 1.5 10"^ < £1 < 2.0 lO^''. We do 
not know any thermal model that predicts a compression 
of the shocked ambient medium behind the shock front in 
such a thin layer, as noted by Vink & Laming (2003) for 
Cas A. Hence, it is likely that the X-ray emission arising 
from this filamentary region is nonthermal in origin. 

For the shocked ejecta, the equality between T2 (cf. 
Eq.El) and T2 (cf. TableEj) implies that £2 = 2.1 10"^ i^. 
Thus, these ejecta fill a very small fraction of the volume 
of the region Al. An upper limit (lower) on £2 is obtained 
by taking the maximum (minimum) values of f2 and the 
minimum (maximum) values of T2: 1.4 10^^ = £2 — 
£2 < £^ = 1.1 10-"* ^2 at a 90% confidence level. These 
values must be compared to £2,chcvaiicr — 0.08 calculated 
from the self-similar model of Chevalier (1982) of n = 9 
(power-law index of the initial outer density profile in the 
ejecta), and s = (power-law index of the initial ambient 
medium density profile). In Appendix IB. 21 Table IBTI indi- 
cates that other values of n and s do not change £2,Chcvaiicr 
a lot. With C - 1/3, we have £2 = 2.3 10"^ (£^ = 1.6 10"^ 
£^ = 1.2 10^^). It is hard to imagine how the ejecta could 
be so clumpy. With such small volumes, the density would 
be very high {he ^ 80 cm~^) and would induce short cool- 
ing timescales (^ 150 years) leading to observable optical 
knots in the southeastern region of Kepler's SNR which 
does not seem to be the case (Bandiera & van den Bergh 
1991; Douvion et al. 2001). Consequently, this model with 
two NEI components does not lead to a consistent astro- 
physical model for the spectrum of the region Al. 

6.2.2. Nonthermal model 

We consider now models for which the emission from the 
shocked ambient medium is nonthermal. The ejecta com- 
ponent of this model is the same as in the previous model. 
As for the shocked ambient medium modelling, rather 



than choosing a simple power-law, we used a synchrotron 
radiation cutolf model (called SRCUT in XSPEC, see 
Reynolds & Keohane 1999) which allows us to describe 
the synchrotron spectrum including the radio data. This 
nonthermal model embraces the following parameters: the 
1 GHz radio flux density 5*1 ghz, the radio spectral index 
a and the roUoff frequency i^r (i-e. approximatively the 
frequency at which the flux has dropped by a factor of 10 
from a straight power-law). 

To estimate the local radio flux 5*1 qhz at 1 GHz in the 
region Al, we used the local radio flux at 1.35 GHz in the 
region Al which is estimated to ~ 0.13 Jy from the radio 
map presented in Dickel et al. (1988). To convert the flux 
at 1.35 GHz to the value at 1 GHz, we normalized it by the 
ratio of the total radio flux at 1 GHz of 19 Jy (Green 2001) 
with the total radio flux at 1.35 GHz of 15 Jy (DeLaney 
et al. 2002), giving S*! ghz = 0.16 Jy. 

Then, we froze the 1 GHz radio flux at 0.16 Jy in 
XSPEC whereas i^r and a were freely fitted. This results 
in a better fit (x^ = 627 for 602 dof) than the fit ob- 
tained with the previous model with 2 NEI (x^ = 701 
for 601 dof). Figure |S1 (right panel) shows the best-fit ob- 
tained with the SRCUT+NEI model and illustrates the 
decomposition of the different components. The best-fit 
parameters are shown in Table [3 The fitted value of the 
radio spectral index a ~ 0.66 ± 0.01 can be compared 
to the global spectral index of 0.64 ± 0.02 (Matsui et al. 
1984) found for the entire SNR and to an estimate of the 
local value of about 0.70 ± 0.02 taken from the spectral 
index radio map of Kepler's SNR made by DeLaney et 
al. (2002). Our value of a is between these values and 
particularly, the derived upper limit of our spectral in- 
dex a — 0.67 is flatter than DeLaney et al.'s value. We 
check that flxing the radio flux to 0.16 and the spectral 
index a to 0.70 makes the flt worse (x^ = 1315). The dif- 
ference between our best flt spectral index and the local 
(spatially and spectrally) radio index can be understood 
because we measure the mean slope between the radio and 
the X-rays. We note that such a curvature of the spectrum 
(radio index steeper than mean index) is predicted by the 
non-linear shock acceleration theory (Reynolds & Ellison 
1992). 

Hence, it is likely that the emission coming from the 
forward shock is nonthermal in origin. However, like in the 
model with two NEI components, the model representing 
the shocked ejecta is still not consistent. 

The spectrum of the forward shock in Kepler is similar 
to that observed by Chandra in Tycho (Hwang et al. 2002) 
and Gas A (Vink & Laming 2003). Using a different argu- 
ments, it was also argued in these papers that a significant 
part of the spectrum behind the forward shock was consis- 
tent with a nonthermal interpretation. For the future, the 
combination of the Chandra and XMM-Newton data will 
provide a deeper understanding and stronger constraints 
on the regions dominated by the forward shock emission. 

The maximum electron energy Emax can be deduced 
from the roUoff frequency Vr and the magnetic field B 
(see Reynolds & Keohane 1999). Using Table we get 



Smax (to|g) = 65^7^ TeV. A magnetic field oi B = 

14 /^G (Matsui et al. 1984) gives E^^^ = 55l " TeV. This 
local value is consistent with the 50 TeV maximum elec- 
tron energy found in Reynolds & Keohane (1999) for the 
fuU SNR. 

6.3. Three components models 

In the previous nonthermal model fSect. IH!2.2() . the ther- 
mal emission from the shocked ambient medium is not 
modelled but it must be there. It is possible to estimate 
its density by adding a third component. 

We add another NEI model to this NEI+SRGUT 
model by fixing all the elemental abundances to their solar 
values. Moreover, the ionization age ri and the normaliza- 
tion Normi of the shocked ambient medium can be com- 
puted respectively from (|B.2|I and from the normalization 
formula ljB.3|l as a function of the mean H density fin. 
Here, we assume that ^ = 1/3 and the parameter ei is 
fixed to ei, Chevalier — 0.3, value chosen from the Chevalier 
(1982) self-similar profiles (n = 9, s = 0). Once the 
shocked gas density fin is fixed, Ti{fiu) and Normi(fiYi) 
are frozen and the data are fitted. We have included in 
XSPEC an upper limit on the temperature Ti given by the 
relation = mnv? with = ~ 4550km/s (cf. 

Appendix IB|| and then we set Ti < Tg — 25 keV. It was 
found that the third component does not improve the fit. 
Then the plot of the curve as a function of nu gives an 
upper limit on njj, i.e. on the shocked ambient medium 
density. We find fi-a < 0.15 cm""^ at a 90% confidence 
level which indicates that the forward shock propagates 
in a very tenuous ambient medium in the southeast. 

7. Conclusion 

The XMM-Newton observation has yielded several new 
results on Kepler's SNR. Observational facts and their as- 
trophysical interpretations are summarized here: 

— The Si K line emission extends to a higher radius than 
the Fe L line emission, notably in the southern part of 
the remnant. It indicates that there is no inversion of 
the Si and Fe layers. 

— The Fe K line emission peaks at a smaller radius than 
the Fe L line emission in the north. It is a clue for a 
higher temperature inwards the remnant. 

— The 4-6 keV and 8-10 keV images show that there are 
few variations in the hardness of the continuum over 
the remnant, except in a few features. 

— The asymmetry in the Fe K emission-line is not associ- 
ated with any asymmetry in the Fe K equivalent width 
map. Hence, the abundance variations do not cause the 
north-south brightness asymmetry. The strong emis- 
sion in the north may be due to overdensities in the 
ambient medium. As for the southeastern region of the 
remnant, it has a very low equivalent width. The Si K 
maps lead to the same conclusions. 



— ihe spectrum ot the southeastern rim indicates that 
the emission of the shocked ambient medium is hkely 
to be mostly nonthermal there. If this is true, the H 
density of the shocked ambient medium in that region 
must be lower than 0.15 cm^'^. 
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Appendix A: Optimization methods for images 

A.l. Optimization of the energy bands width 

The width of the energy bands has been optimized by the 
following method : each line and the surrounding contin- 
uum are fitted by the following model (in cts/keV): 

^{E) ^ Pp{E) -f i?exp (- ^ + R{E) (A.l) 
V ' 

G(E) 

where Pp{E) = A exp{-PE) or P/siE) = AE''^ according 
to the shape of the continuum. R{E) is a reference back- 
ground (Lumb et al. 2002) which has been renormalized in 
the 10-12 keV band for MOS and 12-14 keV band for pn. 
The coefficients A, /3, B,Eo,a of the model are evaluated 
by minimizing the Cash statistic (Cash 1979). Thus, the 
width [El , E2] of each band is determined by optimizing 
the signal to noise ratio: 



S/N{Ei,E2)o, I G{E)dE/[l -^{E)dE\ 
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{E)dE I (A.2) 



A.2. Continuum subtraction 

We can take advantage of the above method to clean line- 
images from the continuum contribution. 

Our goal is to obtain an image of the continuum 
{Ic) which lies under the emission-line in the energy 
band [i?i,i?2] that will be subtracted from the initial 
emission-line image. For that, we use two continuum bands 
and [E^,E2] located on either side of the fine 
band. These two bands are chosen so that the continuum 
is not contamined by nearby lines. Due to the poor statis- 
tics, it is not possible to have directly the shape of the 
continuum in each pixel (i, j). So, we first assume that the 
shape of the continuum remains constant over the SNR. 

Then, the underlying continuum image Ic in the band 
[i?i,i?2] will be a linear combination of the background- 
subtracted continuum images /~ and in the bands 
[E^ ^E^] and respectively. Formally, 



where the normalisation coethcients (a~,aj) and the 
weight A G [0,1] are determined from the spectrum in- 
tegrated over the entire SNR as described below. 

We compute once and for all the normalisation coef- 
ficients by a~ = Mpg/V^ and = Mpg/V^ where V~ 
(I?+) is the measured background-subtracted continuum 
counts in the band [£'f,i?;7] {[^i^^t])- The number of 
counts M[3 in the modelled continuum of slope [3 is defined 
by Pf}{E)dE in the energy band [Si, S2] (A^^ and7V+ 
in the bands [E^,E2] and [E^.E^]^ respectively). The 
value /3o is determined from the adjustement of the entire 
remnant spectrum in the band [i?]^,i?^] as described in 
Sect. lA.ll If the shape of the continuum is the same over 
the entire SNR then, whatever A, the underlying contin- 
uum is correctly represented by Ic defined in Eq. I|A.3|I . 

However, it is possible that the shape of the contin- 
uum changes over the SNR. Then, our goal is to choose 
the weight A to be the least sensitive to spectral varia- 
tions, i.e. to variations of /3 within the model (|A.1|I . In 
each pixel (i, j), we suppose that the modelled continuum 
has a slope (= /3 hereafter). Then, in each pixel, 

the proportions between Ic, /~, I^ should follow those 
between Af/j, Af^,Af^. But, from Eq. HA.3|) . Ic will follow: 



Cfl = A I 



AA,7 + (1-A) a+AA+ 



(A.4) 



By construction, we have CjSg = Npg V A, if the spectral 
model Pp is a good representation of the data. To ensure 
that this equality remains true for small spectral varia- 
tions, we choose the weight A such that: 

{d{Cp-Afp)/dP)p^p^^Q (A.5) 

However, the solution A of Eq. I) A. 5(1 may give too much 
weight to the continuum map which has the lower statis- 
tics so that statistical fluctuations are amplified. To avoid 
this effect, we take the value of A given by Eq. (TOll pro- 
vided that the variance of Eq. (|A.3|I is lower than 10% of 
the variance of the initial emission-line. This is equivalent 
to: 



Vc = [Aa- 



[(l-A)«+]2Ar+ < 0.1 Mine (A.6) 



Ic{i,j) = Xa^ IciiJ) + (1- A)a+/+(z,j) 



(A.3) 



where N" {N^) is the total number of counts {i.e. 
with the background) in the continuum band [E^ , E2] 
{[E^ , E2]) and iVjine the total number of counts in the line 
in [Ei,E2]. If the procedure gives higher variance, then we 
choose the closest A compatible with the constraint on the 
variance {i.e. we give higher weight to the continuum band 
with better statistics). 

A.3. Equivalent width image 

The interpolation based method can also be used to con- 
struct EQW images with some variants. 

To build an EQW image, we have divided the image 
in the line (background and continuum subtracted) by the 
underlying continuum Ic (background subtracted) at Eq. 

Because of such an operation, the statistical fluctua- 
tions can be strongly amplified. So, we must insist on the 



statistical aspect rather than on the spectral aspect, in 
that case, the parameter A of the Ic image used for divid- 
ing is calculated so that the relative error \/V^/Afi3(, (see 
Sect. lA. 2)1 is at most of the same order as the relative error 
on the initial emission- line spectrum l/^/Nunc 

Appendix B: Mean ionization age computation 

Here, we present the method used for computing a mean 
ionization age in the aim to check the value of the ioniza- 
tion age obtained by the spectral fit. 

B.l. General method 

We can write the mean ionization age f as a function of 
the mean electron density fie which is determined from 
the observation, the age of the SNR tg and a parameter ^ 
depending on the spatial distribution, so that 



snocKea amoient meaium 

S £1, Chevalier 



snocKea ejecta 

^2 £2, Chevalier 



(B.l) 



We first calculate the mean electron density as a 
function of the mean H density nu by : 



= C^ZiAbiAbQ^i)n-ii 
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Table B.l. Normalized emission-measure averaged ion- 
ization parameters (^1 and ^2) and filling factors 
(ei, Chevalier and e2,chovaiior) in the self-similar model of 
Chevalier (1982) for different values of n and s. n is the 
power-law index of the initial outer density profile in the 
ejecta and s the power-law index of the initial ambient 
medium density profile. 5 = corresponds to an homo- 
geneous ambient medium, s = 2 to an ambient medium 
modified by a wind and s = — 1 simulates a dense shell in 
the ambient medium. The parameters ^1 and ^2 are calcu- 
lated using Eqs. (|B.4(I and ljB.5|l . It is not possible to give 
the value of ^1 and ^2 for s = 2 because in that case both 
ionization age and density profiles increase to infinity as 
they approach the contact discontinuity. 



where Zi is the number of electrons lost by the atom i 
(fully or partly ionized), Abi and AbQ^i respectively the 
abundance with respect to solar used by XSPEC and the 
solar abundance^ of the element i, i={H,He,C,N,0,Ne, 
Mg,Si,S,Ca,Fe,Ni}. Then, from the normalization formula 
given in XSPEC, we can deduce the mean H density: 



47r 



Norm 



10- 



V z,Ab,Ab 



1/2 



(B.3) 



where D is the distance to the source (in cm), Norm the 
normalization value given by XSPEC and V the emission 
volume which projects onto the local extraction region (we 
shall explain below how to estimate the local volume V). 
Hence, fie is calculated with ljB.2(l using l|B.3|l . 

If the spatial structure is known (from a model) , it is 
possible to evaluate ^. The ionization age should be esti- 
mated by averaging over emission measure (see Kaastra 
& Jansen 1993) 



T(r) rie(r) r dr] / [ ne{r) r dr 



(B.4) 



and the mean electron density from the emission measure 

1/2 



[ / ne{rY r^ dr] / [ r^ dr] 



(B.5) 



so that ^ = T / {fie ts). For reference, f and fie can 
be computed using a self-similar hydrodynamical model 
(Chevaher 1982). In Eq. l|R4jl . T(r) = //^'^^^ ne{t)dt where 
to{r) is the time when a mass element located at the po- 
sition r at tg was shock-heated by the blast wave. In Eqs. 

^ For solar abundances, ZiAbiAbQA — 1.21. 



(IB.4|I and (|B.5|I . the integration is made over the volume 
of the considered shell (shocked ejecta or shocked ambient 
medium) . Table lBTl gives some values of ^ for different val- 
ues of the power-law index n of initial outer density profile 
in the ejecta and the power-law index s of the initial outer 
density profile for the ambient medium. 

B.2. Volume estimation 

Here, we intend to describe how we estimate the emission 
volume V which projects onto the local extraction region. 
This estimation depends on the local geometry of the SNR 
but can be treated with a global model depending on the 
situation encountered. 

Two cases should be considered: 

1, The first one corresponds to the situation in which 
the projected region is contained between the forward 
shock and the reverse shock radii (respectively and 
rr) of the projected shell of the entire SNR. 

2, The second one corresponds to the situation in which 
the shocked ejecta and the shocked ambient medium 
partially fill the volume V which projects onto the local 
extraction region. 

Case 1: The volume V that we integrate along the line-of- 
sight in the local region is always a cylinder full of shocked 
ejecta and shocked ambient medium whatever is the posi- 
tion of the reverse shock. So the fact that we use local or 
global parameters to estimate the volume which projects 
onto the local extraction region does not matter. 



ihc volume V can be computed as a lunction or the 
overall volume of the SNR considered as a filled sphere 
whose radius rg — D Og is given by the position of the for- 
ward shock ^. We define a so that V = a ^ir {D 6sY- 
Afterwards, we correct 1^ by a filling factor e < 1 
which takes into account the volume fraction filled by the 
shocked gas (ambient or ejecta). 

In this configuration, the mean ionization age using 
Eqs. ||RT|) . ||R2)| and (|R3|) becomes: 



r = 9.2 10" ( — , , „ 



Norm ZiAbiAbQ^i 



-3/2 / 2) \ "1/2 



1 kpc 



1/2 



(e a) 



— scm"^ (B.6) 



Case 2: In that case, Eq. HB.6|I is still valid but we ap- 
ply a filling factor e which takes into account the average 
properties of the SNR. 

Table IrTI gives some values of filling factors calculated 
in the framework of the self-similar model of Chevalier 
(1982) averaged over the remnant. Let Ts be the radius 
of the forward shock, r^. the radius of the reverse shock 
radius and Tc the radius of the contact discontinuity. The 
filling factor of the shocked ambient medium is given by 
the geometry of the shell and is equal to £i,chcvaiicr — 
(r^ — rD/r^. On the other side, the filling factor of the 
shocked ejecta is equal to e2,Chevaiier = {^1 - rf)/r^. 

B.3. Recently shocked gas 

The thermal fits (cf. Sect. lt).2TT)) result in very small ion- 
ization age r. This could happen if the X-ray gas was very 
recently shocked. 

We assume that the speed of the shock Vg is constant 
on the time scale ^ tg where tg is the SNR age. Let no 
be the mean ambient density and the compression ratio 
(r2 > 4) . The amount of shocked ambient medium matter 
in the region Al is U,nQeV where V is the volume that we 
integrate along the line-of-sight in the considered region 
and e a fiUing factor (see Sect. IB.2p . It is also S Vg ^tg uq 
where S is the surface encountered by the shock in the 
region Al. By taking an upper limit S'max on the surface 
shocked by the blast wave, we obtain a lower limit ^rnin on 
^. With S'max = 4 TAi rg where tai = D 9ai is the radius 
of the projected region Al and Vg the shock radius, we 
obtain: 



Un = 4.9710-^ 



Ikpcj ^1'7 ^ 1" 



X 



1 yr y \ 1 km/s 



a 



n e 
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® For Kepler's SNR, we set Os = 120" from the 4-6 keV 
continuum image by choosing the same coordinates center as 
in Hughes (1999). This angular radius corresponds to Vs = 
2.8 ±0.8 pc. 



i his equation shows that ^ is geometrically related to the 
filling factor e (cf. Sect. IB.2)l ''. A small value of ^ implies 
a very small filling factor especially as we took a rough 
upper limit on the surface S. 
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